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ABSTRACT 


An  experimental  apparatus  has  been  constructed  and 
tested  to  examine  the  convective  heat  transfer  in  straight 
and  curved  ducts  of  rectangular  cross-section.  The  channel 
can  be  configured  with  both  walls  independently  heated  at  a 
constant  heat  flux,  or  one  wall  at  a  constant  heat  flux  and 
the  opposite  wall  adiabatic.  Local  and  average  Nusselt 
numbers  can  be  calculated,  and  used  to  evaluate  the  effect 
of  curvature  on  heat  transfer. 

Experiments  were  conducted  at  steady  state  for  turbulent 
flow  with  one  wall  heated  at  a  constant  heat  flux  and  the 
opposite  wall  adiabatic.  The  beat  transfer  characteristics 
of  the  straight  and  curved  sections,  on  both  the  inner  and 
outer  walls,  were  compared.  The  heat  transfer  rate  of  the 
concave  curved  wall  proved  to  be  enhanced  over  both  the 
convex  curved  wall  and  both  straight  sections. 
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I. 


INTRODUCTION 


A.  TAYLOR-GORTLER  VORTICES 

Since  the  early  part  of  the  twentieth  century, 
considerable  research  has  shown  that  the  fully  developed 
laminar  flow  along  a  concave  wall  does  not  remain  two- 
dimensional  [Refs.  1,  2,  3].  The  flow  instead  forms  a 
system  of  spiral  vortices,  of  counter  rotating  pairs,  whose 
axes  are  aligned  in  the  principle  flow  direction.  This 
phenomenon  is  the  result  of  the  variations  in  the 
centrifugal  forces  acting  on  the  fluid  particles,  and  is 
known  as  Taylor-Gortler  vortices.  Figure  1  illustrates  the 
type  of  fluid  motion  just  described. 

In  a  channel  that  is  curved  in  the  streamwise  direction, 
those  fluid  particles  located  in  the  center  of  the  flow 
cross-section  are  subject  to  higher  centrifugal  forces  than 
those  fluid  particles  traveling  along  the  channel's  outer 
boundary  wall.  As  a  result,  the  fluid  in  the  center  of  the 
channel  moves  outwardly  toward  the  concave  boundary.  As 
the  process  continues,  the  fluid  particles  near  the  boundary 
wall,  move  in  a  spanwise  direction,  and  finally  radially 
inward  replacing  the  outwardly  moving  particles.  These 
particles  then  come  under  the  same  centrifugal  force  and  the 
process  repeats  itself.  The  resulting  cyclic  motion  causes 
the  formation  of  the  counter  rotating  Taylor-Gortler 


vortices,  considered  to  be  primarily  a  laminar  flow 
phenomenon. 


It  has  been  observed  [Ref.  4]  that  the  heat  transfer 
rate  from  flow  along  a  concave  curved  wall  is  greater  than 
that  for  flow  along  a  straight  wall,  a  phenomenon  attributed 
to  the  additional  mixing  provided  by  the  secondary  motion  of 
the  Taylor-Gortler  vortices. 

There  are  many  possible  applications  that  could  result 
from  a  more  thorough  understanding  of  Taylor-Gortler 
vortices  and  their  effect  on  heat  transfer  and  fluid  flow 
characteristics.  Two  such  applications  could  lead  to 
improved  heat  exchanger  designs  and  improved  turbine  blade 
cooling  [Refs.  5,  6,  7]. 

B.  HISTORY 

The  stability  of  an  inviscid  fluid  flowing  past  a  curved 
boundary,  was  first  considered  by  Lord  Rayleigh  in  1916 
[Ref.  8].  By  assuming  that  the  fluid  was  non-viscous,  he 
determined  that  for  the  motion  to  remain  stable,  its 
circulation  must  increase  with  increasing  radius.  G.  I. 
Taylor,  in  1923,  [Refs.  1,  9],  continued  this  study  with  an 
extensive  analytical  and  experimental  study  of  viscous 
fluids.  His  investigations  focused  on  the  flow  between  two 
cylinders,  in  which  the  inner  cylinder  rotated  while  the 
outer  cylinder  remained  stationary.  Taylor  ascertained  that 
such  flows  become  unstable  when  the  value  of  a  dimensionless 
parameter  exceeded  a  critical  value  of  41.3. 


The  parameter,  known  as  the  Taylor  number  is  defined  as: 

Ta  =  Re 

where  'd'  is  the  width  of  the  gap,  assumed  small  when 
compared  to  'Ri',  the  radius  of  the  inner  cylinder,  and  'Re' 
is  the  Reynolds  number  based  on  the  peripheral  velocity  of 
the  inner  cylinder.  Taylor  determined  that  for  those  cases 
in  which  the  value  of  the  Taylor  number  exceeded  the 
critical  value,  a  secondary  motion  developed  and  the  Taylor 
vortices  formed.  Figure  2  illustrates  this  fluid  motion. 

Instability  of  a  similar  nature  is  also  observed  when  a 
viscous  fluid  flows  in  a  curved  channel  due  to  a  pressure 
gradient  acting  along  the  channel  wall.  This  problem  was 
first  considered  analytically  by  W.  R.  Dean  [Ref.  10]  in 
1928,  for  a  channel  formed  by  two  concentric  cylinders, 
where  the  radius  of  the  inner  cylinder  was  large  in 
comparison  to  the  spacing  between  the  inner  and  outer 
cylinder  walls.  Dean  concluded  that  there  would  be  an 
initiation  of  flow  instability,  and  the  formation  of 
vortices,  when  a  dimensionless  parameter,  the  Dean  number, 
exceeded  a  value  of  36.  The  Dean  number  is  defined  as: 

De  =  Re 

where  'd'  represents  the  channel  half-width,  'Ri '  is  the 
inner  cylinder  radius,  and  'Re'  is  the  Reynolds  number  based 


on  the  mean  velocity  of  the  undisturbed  flow  and  the  channel 
half-width  'd'.  The  analytical  work  of  Dean  was  later 
verified  by  W.  H.  Reid  [Ref.  11],  using  an  approximate 
solution. 

In  1940,  H.  Gortler  [Ref.  2]  studied  the  stability  of 
laminar  boundary  layer  profiles  on  curved  walls  under  the 
influence  of  small  disturbances.  He  found  that  these 
disturbances  were  similar  to  the  vortices  studied  by  G.  I. 
Taylor.  Using  approximate  numerical  calculations,  Gortler 
concluded  that  the  disturbances  of  vortices  were  produced 
only  on  the  concave  boundary  walls  and  that  the  overall  flow 
profile  appeared  to  remain  laminar  in  nature.  These  results 
were  verified,  with  an  exact  solution,  by  G.  Hammerlin,  and 
reported  by  H.  Schlichting  [Ref.  12]  in  1955.  A.  M.  0. 

Smith  completed  an  even  more  extensive  numerical  analysis 
that  further  substantiated  these  findings  [Ref.  3].  The 
results  of  these  numerical  solutions  have  recently  been 
demonstrated  with  the  use  of  hot  wire  anemometry,  laser 
doppler  systems,  and  flow  visualizations  techniques,  [Refs. 
13,  14];  and  in  1976,  Y.  Aihara  [Ref.  15]  conducted  a  non¬ 
linear  analysis  of  Gortler  Vortices. 

As  interest  began  to  develop  concerning  the  effects  of 
the  secondary  flows  associated  with  the  Taylor-Gortler 
vortices,  many  studies  were  published  concerning  the 
influences  of  these  vortices  on  heat  transfer  in  the  laminar 


and  turbulent  flow  regimes.  In  1955,  F.  Kreith  [Ref.  16], 
studied  the  influence  of  heat  transfer  with  respect  to  the 
curvature  of  the  boundary  wall  for  fully  turbulent  flows. 

He  concluded  that  the  heat  transfer  from  the  heated  concave 
boundary  wall  was  considerably  higher  than  that  transferred 
from  the  convex  boundary  wall  of  the  same  curvature  and 
under  similar  turbulent  conditions. 

In  1965,  L.  Persen  [Ref.  17],  considered  the  special 
cases  of  very  high  and  very  low  Prandtl  number  fluids  and 
related  the  increase  in  heat  transfer  from  a  curved  wall  to 
the  presence  of  the  Taylor-Gortler  vortices. 

There  has  been  only  a  limited  amount  of  published 
literature  dealing  with  the  flow  and  heat  transfer  in  curved 
channels  of  rectangular  cross-section  and  large  aspect 
ratios,  where  the  aspect  ratio  is  equal  to  the  spanwise 
distance  of  the  channel  divided  by  the  channel  height.  Much 
of  what  has  been  published  involves  the  development  of 
numerical  approximations  and  solutions  for  heat  and  mass 
transfer  in  curved  ducts  of  various  geometries.  K.  Cheng 
and  M.  Akiyama  [Ref.  18]  developed  a  numerical  solution  for 
forced  convection  heat  transfer  with  laminar  flows  in  curved 
channels  of  rectangular  cross-section,  but  only  for  small 
aspect  ratios. 


In  1976,  A.  A.  Shibani  and  M.  N.  Ozisik  [Ref.  19],  using 
matched  asymptotic  expansion  techniques  for  a  wide  range  of 


Prandtl  numbers,  solved  the  heat  transfer  problem  between 
straight  parallel  plates  with  turbulent  flow,  for  the  case 
of  uniform  wall  temperature. 

Other  experimental  and  analytical  studies  that  are  worth 
citing  with  regard  to  this  present  study  are  Y.  Mori  [Ref. 
20],  who  obtained  results  for  hydrodynamically  fully 
developed  flows  with  constant  wall  heat  flux  in  curved 
channels  of  square  cross-section.  W.  M.  Kays  and  E.  Y. 

Leung  [Ref.  21]  reported  solutions  for  turbulent  flow  heat 
transfer  in  a  concentric  circular  tube  annulus  with  a  fully 
developed  velocity  profile  and  constant  heat  rate  per  unit 
length,  for  a  fluid  of  Prandtl  number  0.7.  P.  F.  Brinich 
and  R.  W.  Graham  [Ref.  22]  reported  results  for  turbulent 
flows  in  a  rectangular  curved  channel  with  an  aspect  ratio 
of  6,  for  the  inner  wall  heated,  the  outer  wall  heated,  and 
both  walls  heated. 

In  1973,  at  the  Naval  Postgraduate  School,  R.  J.  Mckee 
[Ref.  23]  confirmed  the  presence  of  Taylor-Gortler  vortices 
in  a  curved  rectangular  channel  of  aspect  ratio  40.  M. 

Durao  [Ref.  24]  modified  R.  J.  Mckee’s  channel  and  modeled 
it  as  infinite  parallel  plates  with  the  outer  wall  heated 
and  inner  wall  adiabatic.  M.  Durao  and  J.  Ballard  [Refs. 

24,  25]  investigated  the  effects  of  Taylor-Gortler  vortices 
on  the  heat  transfer  in  a  straight  and  curved  test  section 
for  laminar  flows.  R.  Holihan,  Jr.  [Ref.  26]  reported 


results  for  laminar  and  transition  flows.  S.  Daughety 
[Ref.  27]  reported  results  for  turbulent  flows  and  J.  Wilson 
[Ref.  28]  reported  results  for  transition  and  turbulent 
flows.  Each  time  Mckee's  apparatus  was  used  it  was  modified 


slightly  using  prior  recommendations  to  improve  on  the 


accuracy  of  the  results. 


II. 


INTENT  OF  STUDY 


The  purpose  of  this  investigation  was  to  build  an 
improved  curved  rectangular  duct  of  large  aspect  ratio  which 
could  be  configured  with  a  variety  of  boundary  conditions, 
and  yield  more  accurate  data  then  previously  attained.  By 
combining  new  ideas  and  experimental  equipment  with  past 
experiences  and  recommendations  of  R.  Mckee  and  others 
previously  mentioned  [Refs.  23,  24,  25,  26,  27,  28],  the 
effect  of  streamwise  curvature  on  the  heat  transfer  rate  in 
a  curved  rectangular  channel  could  be  studied  in  greater 
detail  while  providing  a  smaller  uncertainty  in  the  results. 

Some  of  the  improvements  incorporated  in  the  new  channel 
include:  the  ability  to  heat  one  or  both  walls  over  a  wide 

range  of  constant  heat  fluxes;  the  ability  to  heat  both 
walls  at  different  constant  heat  fluxes;  the  use  of  more 
calibrated  thermocouples  which  will  allow  greater  precision 
in  reading  surface  and  air  stream  temperatures,  combined 
with  a  better  arrangement  which  will  allow  the  calculation 
of  the  local  heat  transfer  coefficient  and  Nusselt  number; 
a  longer  transition  region  at  the  outlet  which  will  allow 
for  more  mixing  and  a  more  accurate  outlet  temperature 
measurement.  The  use  of  an  outer  and  inner  section  allowed 
both  wall  surfaces  to  be  specially  constructed  with  smooth 


transitions  between  wall  and  heater  surfaces  in  both 
straight  and  curved  test  sections,  preventing  unexpected 
tripping  from  laminar  to  turbulent  flow.  Finally,  by 
constructing  the  channel  with  materials  having  low  thermal 
conductivities  and  avoiding  any  heat  sinks  along  the  channel 
walls,  heat  losses  could  be  minimized. 

Taylor-Gortler  vortices  are  known  to  enhance  heat 
transfer  in  curved  channels  for  transition  and  turbulent 
flows.  This  enhancement  in  the  transfer  of  heat  has  been 
attributed  to  the  secondary  flow  velocity  components  of  the 
Taylor-Gortler  vortices.  These  secondary  components 
transport  the  heated  fluid  from  the  outer  concave  wall, 
inward  toward  the  opposite  wall  of  the  channel,  displacing 
the  cooler  fluid  particles  and  causing  them  to  move  toward 
the  heated  concave  wall.  It  was  expected  that  similar 
results  would  be  observed  in  this  study. 

This  investigation  was  conducted  using  a  single  channel 
with  a  rectangular  cross-section  and  constant  aspect  ratio, 
again  defining  aspect  ratio  as  the  spanwise  distance  divided 
by  the  channel  height.  The  channel  incorporates  a  straight 
and  a  curved  test  section.  Each  test  section  consists  of 
opposing  walls  which  can  be  independently  heated.  The 
results  obtained  in  each  test  section  at  approximately  the 
same  flow  rates  were  compared  in  an  effort  to  determine  the 
effects  of  the  Taylor-Gortler  vortices  on  the  heat  transfer 


rate.  Also,  the  results  of  both  the  straight  and  curved 
outer  sections  were  compared  to  the  results  of  Daughety 
[Ref.  27],  and  Wilson  [Ref.  28]  to  compare  the  new  channel 
to  the  old  one. 

The  straight  section  results  of  this  study  were  compared 
to  the  Dittus-Boelter  correlation  using  hydraulic  diameter 
[Ref.  29],  and  the  results  of  Kays  and  Leung  for  turbulent 
and  transition  flow  in  annular  passages  [Ref.  21]. 

The  curved  section  results  of  this  study  were  compared 
with  the  results  of  Brinich  and  Graham  [Ref.  22]  for 
turbulent  flow  in  a  rectangular  curved  channel,  and  with 
Kreith  [Ref.  16]  for  turbulent  flows  in  concave  and  convex 


curved  channels. 


III.  EXPERIMENTAL  SETUP 


A.  DESCRIPTION  OF  THE  APPARATUS 

To  meet  the  objectives  of  this  investigation,  a  channel 
of  rectangular  cross-section  and  constant  aspect  ratio  was 

I 

constructed.  The  rectangular  channel  has  four  electrically 
heated  wall  sections,  composed  of  a  rubber  heater  bonded  to 
a  thin  copper  plate.  All  four  heated  wall  sections  are 
30.48  centimeters  long  and  25.4  centimeters  wide,  with  a 
heated  area  of  774.2  square  centimeters.  The  rectangular 
channel  is  0.635  centimeters  high  and  25.4  centimeters  wide 
resulting  in  an  aspect  ratio  of  40.  The  cross-sectional 
area  of  the  channel  is  16.13  square  centimeters.  The  wetted 
perimeter  is  52.07  centimeters,  and  the  hydraulic  diameter 
is  1.239  centimeters.  Figure  3  shows  a  cross-sectional  view 
presenting  both  the  straight  and  curved  sections  of  the 
channel . 

The  channel,  shown  in  Figure  4,  has  a  straight  entrance 
region  of  76.20  centimeters,  which  insures  hydrodynamically 
fully  developed  flow  prior  to  entering  the  straight  heated 
test  sections.  The  entrance  region  is  followed  by  two 
parallel  30.48  centimeter  straight  heated  test  sections, 
which  oppose  each  other,  one  in  each  wall.  Following  the 
straight  heated  test  sections  is  a  15.32  centimeter  straight 
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Figure  3.  Schematic  Cross-sectional  View  of  Channel. 


exit  region,  that  completes  the  straight  test  section.  The 
straight  test  section  is  immediately  followed  by  the  curved 
test  section  which  subtends  a  180  degree  arc,  with  the  inper 
wall  having  a  radius  of  curvature  of  29.7  centimeters.  The 
curved  test  section  has  a  52.63  centimeter  unheated  convex 
curved  section  prior  to  the  curved  heated  sections.  The  two 
concentric  curved  heated  sections,  inner  one  convex  and 
outer  one  concave,  are  30.48  centimeters  long  and  the  outer 
heated  wall  subtends  an  arc  of  58.8  degrees.  The  curved 
heated  sections  are  followed  by  a  10.16  centimeter  convex 
curved  exit  region  which  completes  the  curved  test  section. 
Finally,  the  channel  is  completed  by  a  92.0  centimeter 
straight  exit  region  which  allows  sufficient  mixing  prior  to 
reading  the  outlet  temperature.  The  entire  channel  was  then 
thermally  insulated  with  a  plastic  foam  material. 

An  entrance  bell  constructed  of  Plexiglas  is  connected 
to  the  inlet  of  the  straight  entrance  section  of  the 
channel.  It  was  designed  and  manufactured  according  to  ASME 
nozzle  standards,  with  an  elliptical  curved  base  on  a  major 
axis  equal  to  25.4  centimeters  and  a  minor  axis  of  2.54 
centimeters.  A  fiberglass  filter  was  attached  to  the 
entrance  nozzle  to  prevent  foreign  matter  from  entering  the 
channel.  An  aluminum  exhaust  nozzle  was  attached  to  the 
exit  of  the  channel  and  directed  the  flow  from  the  channel 
into  a  two  inch  diameter  pvc  pipe.  The  piping  contained  a 


mass  flow  rate  measuring  section  which  could  be  fitted  with 
standard  ASME  concentric  orifices. 

A  total  of  one  hundred  and  eighty-nine  copper-constantan 
30  gauge  thermocouples  are  located  throughout  the  test 
apparatus.  Only  ninety-nine  thermocouples  are  read  at 
anyone  time,  through  seventy-five  data  channels,  to  record 
the  desired  temperatures.  Each  heated  wall  section  contains 
thirty  thermocouples  which  are  read  individually  and  then 
averaged  to  give  the  wall  section  surface  temperature. 

Four  thermocouples  span  the  width  of  the  channel  at  the 
entrance  and  exit  to  give  the  average  entrance  and  exit  air 
temperatures.  Three  groups  of  five  thermocouples  are 
connected  in  parallel  and  inserted  between  each  layer  of 
insulation  at  each  of  the  heated  test  sections.  Each  group 
of  five  thermocouples  only  uses  one  data  channel,  and  these 
extra  thermocouples  comprise  the  difference  between  the 
number  of  thermocouples  read  and  the  data  channels 
available.  Finally,  there  is  one  thermcouple  to  measure  the 
air  temperature  at  the  orifice. 

B.  DETAILS  OF  THE  CONSTRUCTION 

Each  of  the  four  heated  wall  sections  was  composed  of  a 
copper  plate  laminated  to  a  wire  wound  silicone  rubber 
heater.  The  rubber  heater  was  made  by  Watlow  Industries. 
There  are  thirty  copper-constantan  thermocouples  placed 
between  the  copper  plate  and  rubber  heater.  Figure  5  shows 


Thermocouples 


Outlet 

.egend:  numbers  above  location  marks  are  outer  and 

inner  thermocouple  numbers,  numbers  in 
parentheses  below  location  mark  apply  to 
inner  plates  only. 

Figure  5.  Thermocouple  Placement  in  Copper  Plates. 


the  placement  of  the  thermocouples  on  the  back  of  the  copper 
plate.  The  copper  plate  is  30.48  centimeters  long,  26.72 
centimeters  wide,  and  0.155  centimeters  thick,  but  only  the 
middle  25.4  centimeters  of  copper  plate  width  is  heated  and 
centered  in  the  channel.  The  back  of  the  copper  plate  has  a 
pattern  of  slots  milled  0.076  centimeters  deep.  The  pattern 
of  slots  for  the  two  inner  plates  is  different  from  the 
pattern  of  slots  for  the  two  outer  plates.  See  Figures  6, 

7,  8,  and  9.  The  calibrated  thermocouples  were  numbered  and 
expoxied  into  the  slots,  as  shown  in  Figures  10  and  11.  The 
rubber  heater  has  the  same  length  and  width  dimensions  as 
the  copper  plate,  but  only  heats  the  middle  25.4  centimeters 
of  width.  This  results  in  a  heated  wall  area  of  774.2 
square  centimeters.  The  excess  copper  plate  is  sandwiched 
between  the  wall  and  the  side  spacer  thereby  helping  to 
support  the  copper  plate  in  the  channel.  Dow  Corning  3140 
RTV  was  used  to  adhere  the  rubber  heater  to  the  back  of  the 
copper  plate,  and  complete  the  heated  section.  Figure  12 
and  13  represent  a  completed  test  wall. 

The  channel  was  built  on  a  carefully  constructed  wood 
frame  with  steel  rails,  that  also  was  used  to  support  the 
entire  channel.  The  inner  wall  was  composed  of  two  0.343 
centimeter  thick  and  29.85  centimeter  wide  pieces  of 
Plexiglas  laminated  together,  that  had  spaces  milled  out  for 
the  heated  copper  sections.  The  inner  wall  was  constructed 
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Figure  8.  Photograph  of  Both  Inner  Plates  Showing 
Thermocouple  Slot  Pattern. 


Figure  9.  Photograph  of  Both  Outer  Plates  Showing 
Thermocouple  Slot  Pattern. 


by  heating  and  formi;  the  Plexiglas  around  the  wood  frame. 
The  Plexiglas  wall  was  secured  to  the  frame  by  screws  tapped 
into  the  steel  rails.  The  heated  copper  sections  were  then 
fitted  into  the  machined  spaces  and  two  0.635  centimeter 
thick  and  2.223  centimeter  wide  Plexiglas  spacers,  that 
served  as  the  sides  of  the  channel  and  support  for  the 
copper  plates,  were  attached  by  screws  to  the  inner  wall  and 
steel  rails.  The  two  spacers  were  25.4  centimeters  apart 
and  formed  the  width  of  the  channel.  Additionally,  five 
flat  head  screws  spanning  the  width  of  the  channel  at  the 
joints  where  the  copper  and  Plexiglas  met  were  used  to 
secure  the  copper  plates  and  maintain  an  even  joint.  The 
seams  were  then  worked  to  a  smooth  transition  by  applying 
"Hondo"  to  the  joints  and  sanding.  See  Figures  14,  15,  16, 
and  17. 

The  outer  wall  was  constructed  in  a  similar  manner,  but 
using  0.318  centimeter  thick  and  29.85  centimeter  wide 
pieces  of  Lexan.  Lexan  was  used  on  the  outer  wall  because 
of  its  flexibility  and  ease  of  bending  without  heat,  which 
would  be  required  when  attaching  the  outer  wall  to  the  inner 
wall  and  frame.  Also,  the  outer  wall  was  formed  in  an  outer 
wood  frame  so  that  the  transition  joints  of  the  outer  wall 
could  be  worked  smooth  like  the  inner  wall.  See  Figures  18 
and  19.  Both  walls  were  constructed  so  that  no  seams  or 
butt  joints  existed  on  the  inside  of  the  channel,  except 
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where  the  heated  copper  section  joined  the  Plexiglas  and 
Lexan  walls,  and  all  eight  of  those  joints  had  a  smooth 
transition. 

Once  both  walls  were  completed,  the  outer  wall  and 
support  frame  was  positioned  onto  the  inner  wall  and  main 
frame.  The  outer  wall  was  aligned  onto  the  inner  wall  and 
the  outer  wall  support  frame  removed.  The  outer  wall  was 
then  bolted  to  the  inner  wall  and  frame  with  steel  bars 
running  the  length  of  the  channel  used  as  a  spacer  to  insure 
even  compression.  Once  this  was  completed  the  inlet  and 
outlet  nozzles  were  attached  to  the  channel. 

The  sides  of  the  channel  and  all  connections  between  the 
channel,  the  pvc  piping,  and  the  orifice  were  sealed  with 
Macklanburg-Duncan  Silicone  Rubber  Sealant  Caulk  to  ensure 
there  was  no  leakage  of  air  into  the  channel  or  the  piping 
which  could  effect  the  recorded  temperatures  or  mass  flow 
rate.  The  entire  channel  was  thermally  insulated  with  three 
layers  of  1/2  inch  Armstrong  Armaflex  22  Sheet  Insulation, 
a  flexible  foamed  plastic  material.  By  positioning 
thermocouples  between  these  insulation  layers,  the  heat  lost 
through  the  insulation  could  be  computed  at  each  test 
section.  The  heat  loss  through  insulation  surrounding  the 
rest  of  the  channel  due  to  the  heated  air  flowing  through 
the  channel  was  a  small  fraction  of  the  total  power  supplied 
to  the  heater.  The  insulation  was  held  in  place  by  the  use 
of  silver  duct  tape  and  velcro  strips  attached  to  the  frame. 


c. 


INSTRUMENTATION 


The  data  acquisition  system  used  for  this  experiment  was 
a  Hewlett  Packard  3054A  Automatic  Data  Acquisition/Control 
System  consisting  of  a  3456A  Digital  Voltmeter,  a  3497A  Data 
Acquisition/Control  Unit,  and  four  input  panels  capable  of 
reading  any  mix  of  eight  voltages  or  copper-constantan 
thermocouples.  Also  used  in  conjunction  with  this  data 
acquisition  system,  were  a  Hewlett  Packard  9826  Computer  and 
2671A  Printer.  The  program  used  to  reduce  the  data  is 
listed  in  Appendix  A. 

The  copper-constantan  thermocouples  were  calibrated 
against  a  platinum  resistance  thermometer  using  a  ROSEMONT 
Commutating  Bridge  model  920A,  and  a  ROSEMONT  Constant 
Temperature  Bath  model  913A.  All  thermocouples  were 
calibrated  between  20  and  85  degrees  Celsius,  and  a  second 
order  polynomial  was  curve  fitted  to  correct  each 
thermocouple. 

The  mass  flow  rate  in  the  channel  was  measured  with  a 
thin  plate  orifice  measuring  device  which  was  placed  in  the 
2  inch  pvc  piping  between  the  channel  and  the  compressor. 

In  this  study  an  orifice  with  a  diameter  of  2.7318 
centimeters  was  used  to  measure  hydraulic  Reynold's  numbers 
above  5000.  The  orifice  device  was  installed  in  accordance 
with  ASME  Power  Test  Code  [Ref.  30],  with  3/8  inch  Pressure 
taps  placed  on  either  side  of  the  orifice  at  1  and  1/2 
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diameters,  and  one  thermocouple  inserted  in  the  two  inch  pvc 
piping  down  stream  of  the  orifice  to  record  the  temperature 
of  the  air  flowing  through  the  orifice.  Pressure  measure¬ 
ments  were  taken  with  Meriam  vertical  manometers,  one  whose 
fluid  was  water  with  a  0  to  60  inch  range,  another  using 
water  with  a  0  to  30  inch  range,  and  a  third  whose  fluid 
was  mercury,  calibrated  to  read  inches  of  water  with  a  range 
of  0  to  415  inches.  Atmospheric  pressure  was  measured  with 
a  Princo  precision  barometer.  The  pressure  readings  were 
used  to  calculate  the  mass  flow  rate  of  the  working  fluid  as 
outlined  in  fRefs.  30,  31]. 

Power  was  supplied  to  the  heated  test  section  by  using 
two  Lambda  Model  LK345A  FM  60  volt  power  supplies  in  series. 
The  power  supplies  are  rated  at  6.0  amps,  maximum  at  40 
degrees  Celsius.  Since  the  electrical  resistance  of  the 
rubber  heater  is  not  constant,  but  varies  slowly  with 
temperature,  a  precision  resistor  was  calibrated  and 
connected  in  series  with  the  heater  to  allow  the  calculation 
of  the  instantaneous  power  being  supplied.  Both  the  outer 
and  inner  heated  test  sections  have  their  own  circuitry 
and  precision  resistor.  The  value  of  the  precision 
resistor  for  the  outer  test  sections  is  2.0173  ohms,  and 
the  value  for  the  inner  test  sections  is  2.0084  ohms.  See 
Figure  20  for  diagram  of  circuit  used  to  measure 
instantaneous  power. 


The  working  fluid  was  air  at  room  temperature,  which  was 
drawn  through  the  channel  by  an  electrically  driven  Spencer 
Turbo  Compressor,  rated  at  30  horsepower  at  3500  rpm,  and 
550  cubic  feet  per  minute  at  70  degrees  fahrenheit  and  one 
atmosphere.  Photographs  of  the  channel  and  associated  test 
equipment  are  shown  in  Figures  21,  22,  23,  and  24. 
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IV.  EXPERIMENTAL  PROCEDURES 


The  experimental  procedures  followed  were  the  same  for 
the  straight  and  curved  test  sections.  The  procedure  for 
starting  up  the  experimental  apparatus  is  as  follows: 

(1)  turn  on  the  computer  and  data  acquisition  system 

(2)  close  globe  valve  at  channel  exit 

(3)  close  valve  to  manometer  reading  PI 

(4)  set  compressor  blast  gate  to  1/8  open 

(5)  start  turbo  compressor 

(6)  open  intake  valve  on  compressor  1/4  to  1  turn 

(7)  open  globe  valve  at  channel 

(8)  open  manometer  valve  slowly  to  ensure  PI  does  not 
exceed  manometer  scale 

(9)  turn  on  and  adjust  power  supply  to  the  heater 

(10)  use  intake  valve  on  compressor  to  control  mass  flow 
rate  through  the  channel 

It  was  necessary  to  determine  the  approximate  voltage 
setting  required  to  bring  the  heated  wall  to  approximately 
50  degrees  Celsius.  This  heated  wall  temperature  would 
ensure  a  large  temperature  difference  between  the  inlet  and 
outlet  temperatures,  and  the  heated  and  unheated  wall 
temperatures,  while  not  damaging  the  apparatus.  The  larger 
the  temperature  difference  the  smaller  the  uncertainty  in 
the  results.  The  required  power  setting  was  found  by 


monitoring  the  heated  plate  temperature  locally  from  the 
data  acquisition  system  while  the  voltage  to  the  heater  was 
adjusted . 

To  compute  the  instantaneous  power  supplied  to  the 
heater  the  following  relationship  was  used: 


The  precision  resistor  voltage  and  the  heater  voltage  were 
read  by  the  data  acquisition  system  and  the  resistance  of 
the  precision  resistor  was  a  known  constant.  A  diagram  of 
the  circuit  used  was  shown  in  Figure  16. 

The  mass  flow  rate  was  monitored  by  the  two  manometers 
at  the  orifice.  For  each  run,  the  atmospheric  pressure,  the 
pressure  difference  across  the  orifice,  and  the  pressure 
upstream  of  the  orifice  were  entered  manually  into  the 
computer  for  the  mass  flow  rate  calculation.  All  other  data 
was  acquired  by  the  data  acquisition  system  directly  from 
the  experimental  apparatus. 

Once  the  apparatus  was  started  properly,  sufficient  time 
was  required  to  reach  steady  state  prior  to  gathering  data. 
Preliminary  runs  were  performed  to  determine  the  time 
required  for  the  test  apparatus  to  reach  steady  state.  For 
this  experiment  at  least  two  hours  were  allowed  for  a  cold 
test  apparatus  to  reach  steady  state,  and  at  least  one  hour 


if  the  apparatus  was  already  running.  After  waiting  the 
required  time,  data  was  taken  at  ten  minute  intervals  to 
ensure  that  a  steady  state  condition  had  been  achieved.  Th( 
criteria  for  steady  state  was  based  on  three  variables: 

(1)  the  difference  between  inlet  and  outlet  temperature 

(2)  the  heated  test  section  wall  temperature 

(3)  the  unaccounted  heat  loss  in  the  system 

When  the  first  two  of  these  quantities  varied  by  less  than 
two  percent  over  a  ten  minute  interval  and  the  third  was 
less  than  ten  percent  of  the  input  power,  it  was  considered 
that  a  steady  state  condition  had  been  achieved. 

Experiments  were  conducted  at  specific  mass  flow  rates 
of  air  corresponding  to  Reynolds  numbers  between  7000  and 
22000.  The  Reynolds  number  was  based  on  the  hydraulic 
diameter  of  the  channel.  For  each  mass  flow  rate,  tempera¬ 
ture  and  voltage  data  were  taken  automatically  by  the  data 
acquisition  system  and  immediately  reduced  and  printed.  The 
variables  that  were  measured  and  used  in  this  investigation 
were : 

(1)  the  temperature  of  the  air  entering  the  channel 

(2)  the  temperature  of  the  air  leaving  the  channel 

(3)  the  temperature  of  the  heated  wall  in  the  test 

section  or 

(4)  the  temperature  of  the  unheated  wall  in  the  test 
section  (T„.  or 

(5)  the  temperature  between  each  layer  of  insulation  at 
the  test  section  (T,  _) 


(6)  the  temperature  of  the  air  at  the  orifice 

(7)  the  voltage  across  the  precision  resistor  (Vp^) 

(8)  the  voltage  across  the  rubber  heater  (Vj^) 

(9)  the  pressure  upstream  of  the  orifice  (PI) 

(10)  the  difference  in  the  upstream  and  downstream 
pressure  across  the  orifice  (^P) 

(11)  the  atmospheric  pressure 

The  results  obtained  from  the  straight  test  sections 
served  as  a  baseline  for  the  comparison  of  the  curved 
section  results.  An  effort  was  made  to  match  the  mass  flow 
rates  for  each  test  section  as  closely  as  possible. 


V.  PRESENTATION  OF  DATA 


A.  DATA  REDUCTION 

The  constant  heat  flux  surface  for  the  heated  wall  of 
each  test  section  was  provided  by  the  uniform  electrical 
resistivity  of  the  rubber  heater.  The  insulated  unheated 
wall  was  considered  adiabatic  since  the  heat  losses  through 
that  wall  were  negligible.  The  straight  portion  of  the 
channel  upstream  of  the  straight  heated  test  section  was  of 
sufficient  length  to  ensure  that  the  air  flow  was  hydro- 
dynamically  fully  developed  for  the  flow  velocities  of  this 
experiment.  The  straight  portion  of  the  channel  downstream 
of  the  curved  section  proved  to  be  of  sufficient  length  to 
ensure  that  the  air  flow  exiting  the  channel  was  thoroughly 
mixed  and  that  the  measured  outlet  air  temperature 
was  an  average  bulk  temperature,  as  long  as  the  hydraulic 
Reynolds  number  of  the  channel  was  less  than  17000. 

The  analysis  of  the  data  required  several  quantities  to 
be  defined.  The  heat  convected  to  the  air  was  calculated 
using  the  equation: 

®air  “  *"^pair  ^"^out  ”  "^in^ 

where  'C__.  '  is  the  specific  heat  of  the  air  at  constant 

pd  1  k 

pressure  and  'm'  is  the  mass  flow  rate  of  the  air. 
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The  average  heat  transfer  coefficient  between  the  heated 


wall  and  the  flow  of  air  in  the  channel  was  defined  by  the 
equation: 


h  = 


®ai  r 


where  'Q,.  '  is  defined  above,  'A_,r  '  is  the  area  of  the 

heated  copper  plate  in  the  test  section,  and  'AT'  is  the 
difference  between  the  arithmetic  average  heated  wall 
temperature  (T^)  and  the  average  bulk  air  temperature 
(Tbii^)  .  The  average  bulk  air  temperature  is  defined  as  the 
arithmetic  mean  of  the  entrance  and  exit  air  temperatures 

<''in' 

The  average  Nusselt  number  was  then  calculated  using: 


h  Dhd 

K  . 


air 


In  this  equation  is  the  hydraulic  diameter  and 

is  the  thermal  conductivity  of  the  air. 

The  hydraulic  Reynolds  number  was  calculated  for  each 
test  run  as  follows: 


«®hd 


m  D 


hd 

air 


where  again  'm'  and  are  the  mass  flow  rate  and 

hydraulic  diameter  of  the  channel,  'A  '  is  the  cross- 


sectional  area  of  the  channel,  and  is  the  dynamic 

viscosity  of  the  air. 

For  the  curved  sections  the  Dean  number  was  defined  as 


where  is  the  Reynolds  number  based  on  the  channel 

height,  is  the  channel  height,  and  'Ri '  is  the  radius 

of  curvature  of  the  inner  convex  wall. 

A  sketch  of  the  control  volume  and  a  set  of  sample 
calculations  for  one  test  run  of  the  curved  section  are 
given  in  Appendix  B.  A  sample  calculation  for  the 
uncertainty  analysis  is  given  in  Appendix  C. 

B.  RESULTS 

The  data  obtained  from  each  experimental  run  was  reduced 
by  the  computer  program  listed  in  Appendix  A,  utilizing  the 
expressions  described  in  Part  A  above.  The  major  parameters 
resulting  from  this  process  are  shown  in  Tables  I  through 
IV.  Tables  I  and  II  contain  the  straight  sections  results 
and  Tables  III  and  IV  contain  the  results  of  the  curved 
sections.  A  plot  of  the  average  Nusselt  number  versus 
hydraulic  Reynolds  number  is  given  in  Figures  25,  26,  27, 

28,  and  29  for  the  comparison  of  the  results.  "Error"  bands 
have  been  indicated  as  a  result  of  the  uncertainty  analysis. 


TABLE  I 


SUMMARY  OF  STRAIGHT  OUTER  WALL  TEST  RESULTS 


''^hd 

h 

(W/mC) 

AT 

(C) 

^%d 

7709 

108. 50 

41.99 

28.  52 

20.07 

11953 

135.26 

61.57 

25.49 

29.51 

15646 

157, 92 

74.  43 

24.  45 

35.  65 

21455 

217.53 

87.42 

27.  36 

41.95 

Reynolds 

number  and 

Nusselt 

number  are  based 

on 

hydraulic  diameter.  Results  for  Reynolds  number  of 
21455  based  on  orifice  air  temperature  instead  of 
outlet  air  temperature. 
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TABLE  II 


SUMMARY  OF  STRAIGHT  INNER  WALL  TEST  RESULTS 


^®hd 

Q  • 

h 

(W/mC) 

AT 

(C) 

»^%d 

7628 

105.57 

41.52 

28.90 

19.85 

12169 

138. 36 

59.  50 

27.  42 

28.  50 

15788 

159.84 

70.24 

26.37 

33.61 

21312 

204.42 

81.  82 

27.  84 

39.17 

Reynolds  number  and  Nusselt  number  are  based  on 
hydraulic  diameter.  Results  for  Reynolds  number  of 
21312  based  on  orifice  air  temperature  instead  of 
outlet  air  temperature. 


TABLE  III 


SUMMARY  OF  CURVED  CONCAVE  WALL  TEST  RESULTS 


^®hd 

De 

h 

(W/mC ) 

T 

(C) 

7782 

206 

115. 88 

52.  43 

25.  51 

25.  00 

11726 

311 

157.45 

68. 16 

29.91 

32.53 

15813 

419 

158. 17 

83. 13 

22.  34 

39.78 

21246 

563 

231.31 

93.02 

27.20 

44.52 

Reynolds 

number 

and  Nusselt  number 

are  based 

on 

hydraulic  diameter.  Results  for  Reynolds  number  of 
21246  based  on  orifice  air  temperature  instead  of 
outlet  air  temperature. 


I 


! 

i 


TABLE  IV 


SUMMARY  OF  CURVED  CONVEX  WALL  TEST  RESULTS 


^®hd 

De 

h 

(W/mC ) 

AT 

(C) 

%d 

7757 

206 

82.37 

38.33 

24.47 

18.37 

12021 

318 

115. 50 

57.  62 

23.  92 

27.65 

15908 

421 

142.54 

68.57 

24.34 

32.94 

21603 

572 

164. 08 

79.28 

23.43 

38.  11 

Reynolds  number  and  NusseXt  number  are  based  on 
hydraulic  diameter.  Results  for  Reynolds  number  of 
21603  based  on  orifice  air  temperature  instead  of 
outlet  air  temperature. 
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Figure  26.  Present  Concave  vs.  Convex  Section  Results 
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Figure  28.  Present  Concave  vs.  Straight  Outer  Section  Results 


Figure  29.  Present  Convex  vs.  Straight  Inner  Section  Results 


As  would  be  expected  the  results  indicate  an  increase  in 
the  rate  of  heat  transfer  with  increasing  Reynolds  number 
for  both  the  straight  and  curved  test  sections.  In 
addition,  the  heat  transfer  rate  was  highest  in  the  concave 
curved  section  for  each  Reynolds  number  investigated,  and 
lowest  in  the  convex  curved  section. 

In  the  low  end  of  the  turbulent  regime,  at  a  hydraulic 
Reynolds  number  of  7780,  the  rate  of  heat  transfer  was  about 
twenty-two  percent  higher  in  the  concave  curved  section  than 
the  outer  straight  section.  As  the  hydraulic  Reynolds 
number  continued  to  increase,  between  11700  and  15800,  the 
rate  of  heat  transfer  was  approximately  eleven  percent 
higher  between  the  ccncave  curved  section  and  the  outer 
straight  section.  At  a  hydraulic  Reynolds  number  of  21000 
the  difference  in  heat  transfer  rate  decreased  to  about 
seven  percent. 

The  results  obtained  at  a  hydraulic  Reynolds  number  of 
21000  have  a  much  larger  uncertainty,  because  the  orifice 
air  temperature  was  higher  than  the  average  outlet  air 
temperature.  This  signifies  that  the  air  was  not  fully 
mixed  at  the  outlet,  because  no  heat  is  added  to  the  air 
after  the  outlet  air  temperature  is  read,  yet  the  air 
temperature  increased.  For  hydraulic  Reynolds  numbers 
greater  than  about  17000  the  channel  outlet  section  is  not 
long  enough  to  allow  the  air  to  mix  properly. 


The  values  obtained  for  the  hydraulic  Reynolds  number  of 
approximately  21000  are  based  on  the  orifice  air  temperature 
vice  the  outlet  air  temperature.  This  resulted  in  a  Nusselt 
number  higher  than  obtained  by  using  the  the  outlet  air 
temperature,  but  still  lower  than  the  actual  Nusselt  number, 
because  of  the  small  heat  loss  between  the  outlet  nozzle  and 
the  thermocouple  at  the  orifice.  This  heat  loss  cannot  be 
accounted  for  and  the  temperature  used  as  an  outlet  air 
temperature  was  still  lower  than  the  actual  outlet  air 
temperature . 

In  comparing  the  concave  curved  wall  to  the  convex 
curved  wall  the  rate  of  heat  transfer  was  about  thirty-six 
percent  higher  in  the  concave  section  for  a  hydraulic 
Reynolds  number  of  7780.  For  hydraulic  Reynolds  numbers 
between  11700  and  15800  the  rate  of  heat  transfer  was  twenty 
percent  higher  in  the  concave  section,  and  at  a  hydraulic 
Reynolds  number  of  21200  the  rate  of  heat  transfer  was 
seventeen  percent  higher. 


VI.  DISCUSSION  AND  CONCLUSIONS 

J.  Wilson  [Ref.  28],  in  his  study,  reported  that  the 
actual  presence  of  Taylor-Gortler  vortices  were  verified  by 
the  observation  of  the  liquid  crystal  isotherm  distribution 
in  the  curved  test  section.  He  further  reported,  that  at 
low  values  of  hydraulic  Reynolds  number  and  through  most  of 
the  transition  regime,  liquid  crystal  bands  were  oriented  in 
the  streamwise  direction.  As  the  Reynolds  number  increased 
turbulent  forces  began  to  dominate  in  the  curved  section  and 
less  evidence  could  be  detected  for  the  existence  of  the 
vortices . 

The  high  aspect  ratio  of  the  channel  provided  the  basis 
for  the  assumption  that  the  experimental  apparatus  should 
behave  like  infinite  parallel  plates.  Holihan's  data  [Ref. 
26]  tended  to  verify  this  assumption.  His  experimental  data 
in  the  laminar  flow  region,  approached  the  theoretical  limit 
for  average  Nusselt  number  for  parallel  plates  with  one  wall 
heated  at  a  constant  heat  flux  and  the  opposite  wall 
adiabatic  [Ref.  32] .  Based  on  the  assumptions  mentioned 
above,  comparisons  were  made  with  other  analytical  and 
experimental  results  that  were  of  the  same  or  similar 
problems . 

A  comparison  of  the  experimental  results  of  this  study 
and  the  experimental  results  of  Daughety  [Ref.  27]  and 
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Wilson  [Ref.  28]  are  shown  in  Figure  30  for  the  curved 
section,  and  Figure  31  for  the  straight  section.  The 
channel  aspect  ratio  was  40  for  each  of  these  studies  and 
the  experimental  procedures  were  similar,  with  the  exception 
of  a  slight  modification  in  Wilson's  procedures.  Wilson 
[Ref.  28]  made  a  correction  to  determine  the  temperature 
drop  across  the  liquid  crystal  sheet  covering  the  heated 
wall.  By  applying  this  temperature  drop  to  the  measured 
temperature,  the  true  temperature  of  the  heated  wall  would 
be  known.  After  reviewing  Wilson's  data  and  experimental 
procedure  it  was  detected  that  Wilson  was  over  estimating 
the  temperature  drop  across  the  liquid  crystal  by 
approximately  five  to  seven  percent,  and  this  lead  to  an 
over  estimation  of  T  of  ten  to  twelve  percent.  This 
resulted  in  an  over  estimation  of  the  average  heat  transfer 
coefficient  and  the  average  Nusselt  number  of  between  ten 
and  twelve  percent.  When  Wilson's  data  is  lowered  by  ten 
percent,  both  his  data  and  Daughety's  data  fall  within  the 
uncertainty  bands  of  the  present  study  as  shown  in  Figures 
32  and  33.  Differences  were  attributed  to  the  inability  to 
acccount  for  all  of  the  heat  transfer  processes  and  the 
limitations  inherent  in  any  experimental  work. 

Ballard  [Ref.  25]  reported  an  increase  in  the  rate  of 
heat  transfer  between  a  concave  wall  and  a  straight  outer 
wall  of  eleven  percent  for  laminar  flows.  Holihan  [Ref.  26] 
reported  an  increase  of  fifteen  percent  for  laminar  flows 


70 


LEGEND 


BOVUaAV 


Figure  32.  Comparison  of  Present  Data  with  Daughety  and  Corrected  Wilson,  Curved 
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Figure  34.  Comparison  of  Present  Data  with  Dittus-Boelter ,  Kays  and  Leung 
Turbulent  Flow,  Straight  Section. 


[Ref.  29]  for  heat  transfer  in  a  straight  tube  for  constant 
wall  temperature  is  plotted  from: 

Nu  =  0.02Re°’® 

For  a  Prandtl  number  equal  to  0.71.  Additionally,  the 
results  of  Kays  and  Leung  [Ref.  21]  for  heat  transfer  in 
annular  passages  are  shown  for  r*  =  1.0,  which  equates  to 
parallel  plates.  When  examined  these  results  compare 
favorably  with  the  present  data. 

The  curved  section  data  is  plotted  in  Figure  35,  and  is 
compared  to  the  experimental  work  of  Brinich  and  Graham 
[Ref.  22].  The  accuracy  of  the  data  points  used  from 
Brinich  and  Graham  are  subject  to  some  error,  in  that  the 
actual  values  are  not  given  in  their  study  and  had  to  be 
taken  from  a  plot  of  Stanton  number  versus  Reynolds  number. 

After  carefully  examining  the  past  and  present  data  the 
trends  and  results  are  consistent.  Heat  transfer  rates 
increase  with  increasing  Reynolds  numbers.  Also,  the 
concave  curved  section  enhances  the  rate  of  heat  transfer 
over  the  convex  curved  section  and  the  straight  test 
sections  for  laminar,  transition,  and  turbulent  flows.  The 
presence  of  Taylor-Gortler  vortices  in  the  laminar  and 
transition  regime  seems  to  be  the  major  contributing  factor 
to  this  enhancement.  Whether  the  Taylor-Gortler  vortices 
continue  to  exist  in  the  turbulent  flow  regime  remains  to  be 


HYDRAULIC  REYNOLDS  NUMBER 

Figure  35.  Comparison  of  Present  Data  with  Brinich  and  Graham,  and  Kreith,  for 
Turbulent  Flow,  Curved  Section. 


VII.  RECOMMENDATIONS 


The  amount  of  experimental  work  remaining  in  the  study 
of  heat  transfer  in  curved  channels  is  still  significant. 

The  present  experimental  apparatus  has  been  tested  in  the 
turbulent  regime,  and  the  test  data  proves  that  accurate 
results  can  be  obtained  for  hydraulic  Reynolds  numbers  below 
17000,  slightly  higher  when  heating  the  straight  test 
section.  Additional  experiments  into  the  laminar, 
transition,  and  turbulent  region  should  be  conducted  to 
better  understand  the  heat  transfer  processes  occurring  in 
the  straight  and  curved  sections.  By  remaining  below  a 
hydraulic  Reynolds  number  of  17000  accurate  data  can  be 
obtained  and  significant  conclusions  can  be  drawn  in  the 
comparison  of  the  curved  sections  to  the  straight  section, 
and  the  concave  curved  wall  to  the  convex  curved  wall. 

To  improve  the  accuracy  of  the  data  and  results  some 
slight  modifications  to  the  present  channel  should  be  done. 
First,  the  straight  outlet  section  of  the  channel  should  be 
lengthened  to  ensure  proper  mixing  of  the  heated  outlet  air. 
Second,  some  sort  of  honey  combing,  screening,  or  other  type 
of  device  should  be  placed  in  the  outlet  section  to  enhance 
the  mixing  of  the  heated  outlet  air.  Finally,  extra 
thermocouples  should  be  placed  at  the  outlet  and  read  more 


then  once  during  a  data  run  to  decrease  the  uncertainty  of 
the  outlet  air  temperature  and  the  uncertainty  of  the 

difference  between  the  inlet  and  outlet  air  temperature.  By 
applying  these  recommendations  the  uncertainty  of  the 
Nusselt  number  can  be  reduced  to  four  percent  or  less. 

Further  experiments  at  hydraulic  Reynolds  numbers 
greater  than  17000  also  need  to  be  examined.  This  can  be 
accomplished  if  the  above  recommendations  are  taken  along 
with  one  other.  The  channel  walls  need  outer  support  to 
prevent  the  center  of  the  channel  from  collapsing  at  the 
high  hydraulic  Reynolds  numbers.  Slight  collapsing  was 
noticed  at  hydraulic  Reynolds  numbers  greater  than  20000. 
This  collapsing  was  caused  by  the  laminating  of  two  0.3175 
centimeter  thick  walls  which  were  not  as  rigid  as  the  single 
0.635  centimeter  thick  walls  used  previously.  The  actual 
hydraulic  Reynolds  numbers  that  can  be  attained  and  still 
yield  accurate  results  will  be  limited  by  the  power  of  the 
compressor  and  the  rigidity  of  the  channel's  walls. 

Finally,  once  a  strong  data  base  can  be  obtained  for 
each  test  section  further  experimentation  can  be  conducted 
with  different  boundary  conditions.  For  example,  both 
curved  walls  heated  at  the  same  constant  heat  flux.  Once 
this  can  be  done  a  true  comparison  of  the  heat  transfer  rate 
of  the  total  curved  test  section  can  be  compared  to  the 
total  straight  test  section,  i.e.,  both  curved  walls  versus 
both  straight  walls. 


APPENDIX  A:  COMPUTER  PROGRAM 


1C00  ! 

1010  ! 

1020  ! 

1030! 

1040  ! 

1050  ! 

1060  I 
1 070  ! 

1080  I 
1090! 

1100  COM  /Co/  D(7) ,Aa(7G.2) 

1110  DIM  Eiiif(79),T(79) 

1120! 

1130  !  CORRELATION  FACTORS  FOR  EMF  TO  DEGREES  CELSIUS  FOR  CU-CO  T/C. 

1140! 

1150  DATA  0.10086091,25727.94369,-767345.8295.78025595.81 
1160  DATA  -9247486589, 6. 97688E+11, -2. 66192E+13. 3. 94078E+14 
1170  READ  D(*) 

1180  PRINTER  IS  701 
1190  BEEP 

1200  INPUT  ‘ENTER  MONTH.  DAY,  AND  TIME  (MM:DD:HH:MM:SS) * . TimeS 
1210  OUTPUT  709;'TD‘;Txme$ 

1220  BEEP 

1230  INPUT  ‘ENTER  INPUT  MODE  ( 0-3054A. 1 -FILE) ‘ , In 
1240  IF  lii-l  THEN 
1250  BEEP 

1260  INPUT  "ENTER  NAME  OF  EXISTING  DATA  FILE* ,01df i leS 

1270  PRINT  USING  ■20X, “These  Results  Are  From  Data  File:  “.10A‘:01dfile$ 

1280  PRINT 
1290  PRINT 
1300  ELSE 
1310  BEEP 

1320  INPUT  ‘ENTER  NAME  OF  NEW  DATA  FILE?’ .Newf i leS 

1330  PRINT  USING  •20X,“Data  For  This  Run  Is  Stored  In  Data  Filet  “,10A‘;Nowf 
ileS 

1340  PRINT 
1350  PRINT 
1360  END  IF 
1370! 

1380  !  ENTER  THE  CHANNEL  CONFIGURATION. 

1390! 

1400  PRINT  USING  ‘9X,“THE  FOLLOWING  CHANNEL  CONFIGURATION  HAS  SELECTED:  “* 

1410  PRINT 
1420  BEEP 

1430  INPUT  ‘SELECT  HEATER  TYPE  (0-STRAIGHT,  1-CURVED) ‘ , Itype 

1440  IF  Itype-0  THEN  PRINT  USING  •12X, “Heating  STRAIGHT  Test  Section."* 

1450  IF  Itype-1  THEN  PRINT  USING  •12X,  “Heatinq  CURVED  Test  Section. "“ 

1 460  BEEP 

1470  INPUT  ‘SELECT  HEATER  POSITION  (O-OUTER.  1 -INNER) ‘ . Ipos 
1480  IF  Ipos-O  THEN  PRINT  USING  ‘12X, “Heating  OUTER  Plate. “■ 

1490  IF  Ipos-1  THEN  PRINT  USING  ‘12X, “Heating  INNER  Plate. ‘“ 

1500  PRINT 
15101 

1520  !  ENTERS  FILES  WITH  T/C  CALIBRATION  COEFFICIENTS  BASED  ON  HEATER  TYPE. 

1530! 

1540  IF  Itype-0  THEN  ASSIGN  @Filec  TO  ‘STRA_COEF‘ 

1550  IF  Itype-1  THEN  ASSIGN  @Filec  TO  ‘CURV_C0EF‘ 

1560  FOR  I-O  TO  76 

1570  ENTER  @F i led Aa ( I . 0 ) , Aa(I . 1 ) . Aa( I ,2) 

1580  NEXT  I 
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PROGRAM  FOR  GATHERING  AND  REDUCING  DATA  WHEN  HEATING  ONE  COPPER  PLATE, 
WITH  10  MINUTES  BETWEEN  RUNS.  PROGRAM  IS  MODIFIED  TO  ACCOUNT  FOR 
BROKEN  THERMOCOUPLE  CHANNELS  5.  16,  21.  AND  26  WHICH  EFFECT  OUTER 
CURVED  HEATER,  AND  THERMOCOUPLE  CHANNEL  47  WHICH  EFFECTS  INNER  CURVED 
HEATER. 
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1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

1860 

1870 

1880 

1890! 

1900 

1910! 

1920 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

dd) 

2000 

dd) 

2010 

s)"" 

2020 

s)... 

2030 

2040 

2050 

2060! 

2070 

2080! 

2090 

2100 

2110 

2120 

2130 

2140 


ASSIGN  @Filec  TO  * 

!  INSERT  DATA  STORAGE  DISK. 

BEEP 

INPUT  "CHANGE  DISK  AND  HIT  ENTER", Ok 

IF  In-I  THEN  ASSIGN  SFile  TO  OldfileS 

IF  lm-0  THEN 

CREATE  BOAT  Neufile$,40 

ASSIGN  »File  TO  NewfileS 

END  IF 

ENTER  709; Times 
CLEAR  709 

PRINT  USING  "20X, ""Month,  Day,  and  Time:  "", 14A";TimeS 
PRINT 

!  FOR  DATA  RUN  READ  AND  STORE  ALL  RAM  EMF  VALUES. 

IF  lm-0  THEN 

OUTPUT  709; "AR  AFOO  AL79" 

OUTPUT  722;"F1  R1  T1  Z1  FL1" 

FOR  I-O  TO  79 
OUTPUT  709; "AS" 

ENTER  722;Emf(I) 

NEXT  I 

OUTPUT  eFile;Emf(*) 

ELSE 

ENTER  eFile;Emf(*) 

END  IF 

OUTPUT  709: "TD" 

!  CONVERT  AND  CALIBRATE  ALL  T/C  READINGS  TO  DEGREES  CELSIUS. 

PRINT  USING  *9X,""THE  FOLLOHING  DATA  HERE  RECORDED:""* 

PRINT 

FOR  I-O  TO  79 

IF  1-38  OR  1-39  OR  1-77  OR  1-78  THEN 

IF  (Ipos-0  AND  1-38)  OR  (Ipos-1  AND  1-77)  THEN 

Iadd-0 

IF  Ipos-1  THEN  Iadd-39 

PRINT  USING  "12X,**Hoater  Voltage  (Vh)  -  "".2D. 3D.""  ( V) "" " ;Emf (38+Ia 

PRINT  USING  "12X. ""Resistor  Voltage  (Vr)  -  "".2D. 3D,""  ( V) " " " ;Emf (39+Ia 

IF  Ipos-0  THEN  PRINT  USING  "12X. ""Precision  Resistor  (Rpr)  -  2.01729  (Ohm 


IF  Ipos-1  THEN  PRINT  USING  " 1 2X, ""Precision  Resistor  (Rpr) 

PRINT 
END  IF 
ELSE 

!  IF  T/C  OUT  OF  RANGE  PRINT  A  OUE. 

Aemf -Emf ( I ) 

IF  Aemf<3,E-4  OR  Aerof>3.E-3  THEN 

T(I)-999.99 

ELSE 

CALL  Tvsv(Emf (I) ,Tt) 

IF  1-79  THEN 


2.01729  (Ohm 
2.00839  (Ohm 


2150  T(I)-Tt 

2160  GOTO  2210 

2170  END  IF 

2180  T(I)-FNTcorr(Tt,I) 

2190  END  IF 

2200  END  IF 

2210  NEXT  I 

2220! 

2230  I  PRINT  ALL  T/C  READINGS. 

22401 

2250  PRINT  USING  ■12X, "OUTER  PLATE:"’ 

2260  PRINT  USING  ’MX.  "T/C  Number:  12  3  4  5 

6*" 

2270  PRINT  USING  ■14X,"Temp  (C):  "  ,S(3D.DD,2X) * ; T{ 0) ,T( 1 ) ,T(2) .T(3) .T(4) , T 

(5) 

2280  PRINT  USriG  ’MX. "T/C  Number:  7  8  9  10  11 

12”. 

2290  PRINT  USING  “^X.^’Iemp  (C):  "  ,6(3D.DD,2X) * ; T(6) , T(7 ) , T (8) .T(9) ,T( 1 0) , 

T(11 ) 

2300  PRINT  USING  ■14X,"T/C  Number:  13  14  15  16  17 

18"* 

2310  PRINT  USING  ■14X,"Temp  (C):  "  ,6(3D.DD.2X) * ;T( 12) , T( 13) . T( 1 4) ,T ( 15) .T( 

16) .T(17) 

2320  PRINT  USING  ’^X.^’T/C  Number:  19  20  21  22  23 

24’" 

2330  PRINT  USING  *14X."Temp  (C):  ”  .6(3D.DD,2X) ’ ;T( 18) .T{ 19) .T(20 ) .T(21 ) , T( 

22),T(23) 

2340  PRINT  USING  "14X,"T/C  Number:  25  26  27  28  29 

30’" 

2350  PRINT  USING  “^X.^’Iemp  (C):  " ,6(3D.DD.2X) * ! T(24) , T(25) , T(26) . T(27) , T( 

28),T(29) 

2360  PRINT 
23701 

2380  I  CALCULATE  AVERAGE  TEMPERATURE  OF  OUTER  HEATED  HALL. 

2390! 

2400  Sum“0. 

2410  IF  Itype-0  THEN 
2420  FOR  I-O  TO  29 
2430  Sam”Sum+T(I) 

2440  NEXT  I 
2450  Two~Sum/30 
2460! 

2470  !  CORRECTION  FOR  BROKEN  T/C  CHANNELS  5,  16,  21,  AND  26  IN  OUTER  CURVED  HE 

A  TER. 

2480! 

2490  ELSE 

2500  FOR  I-O  TO  29 

2510  IF  1-5  OR  1-16  OR  1-21  OR  1-26  THEN  2530 
2520  Sum-Sum-*'T(  I ) 

2530  NEXT  I 
2540  T«o-Sum/26 
2550  END  IF 

2560  PRINT  USING  ■12X. "OUTER  INSULATION:’" 

2570  PRINT  USING  •14X,"T/C  Number:  31  (Tinsol)  32  (Tinso2)  33  (Tinso3)’ 

M  ■ 

2580  PRINT  USING  "14X,"Temp  (C):  "  ,3(3D. DD,8X) " : T (30 ) , T(31 ) ,T(32) 

2590  Tjnso1-T(30) 

2600  Tinso2-T(31) 

2610  Tinso3-T(32) 

2620  PRINT 
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2B30  PRINT  USING  ’IZX, "ORIFICE  TEMP  (Torf)  -  ".3D.DD."  (C)"*;T(33) 

2G40  PRINT 

2650  PRINT  USING  •12X, "INLET  TEMPERATURE:*" 

2660  PRINT  USING  ■1AX,*"T/C  Number!  35  36  37  38"" 

2670  PRINT  USING  "MX. "Temp  (C):  **  ,4(3D.DD,2X)  "  ;T(34) ,  T(35) ,  T(36) .  TC37) 

2680  PRINT 
2690! 

2700  I  CALCULATE  AVERAGE  INLET  TEMPERATURE. 

2710! 

2720  Sum'O. 

2730  FOR  I-3A  TO  37 
2740  Sura“Sum+T(I) 

2750  NEXT  I 
2760  Tin”Sui»/4 

2770  PRINT  USING  "12X." INNER  PLATE:"* 

2780  PRINT  USING  ■14X.*"T/C  Number:  41  42  43  44  45 

46*** 

2790  PRINT  USING  ‘^X.^'Iemp  (C):  " .6(3D.DD,2X) * ; T(40 ) , T(41 ) . T(42) ,T(43) , T( 

44) ,T(45) 

2800  PRINT  USING  ■14X,**T/C  Number:  47  48  49  50  51 

52". 

2810  PRINT  USING  ■14X,**Temp  (C):  " .6(3D.DD,2X) * ;T(46) , T(47) , T(48) ,T(49) , T( 

50),T(51) 

2820  PRINT  USING  •14X,"T/C  Number:  53  54  55  56  57 

58"" 

2830  PRINT  USING  "14X,"'Temp  (C):  "* ,6(3D.DD.2X) " ; T (52) , T(53) , T(54) , T (55) , T( 

56),T(57) 

2840  PRINT  USING  •14X,"*T/C  Number!  59  60  61  62  63 

64""" 

2850  PRINT  USING  “14X."Temp  (C):  "".6(3D.DD.2X)*:T(58) .T(59) .T(60) .T(61 ) ,T( 

62),T(63) 

2860  PRINT  USING  ■14X,""T/C  Number!  65  66  67  68  69  " 

70""" 

2870  PRINT  USING  "14X.*"Temp  (C):  "" ,G(3D.DD,2X) " ; T(64) ,T(65) ,T(6G) ,T(67) ,T( 

68),T(69) 

2880  PRINT 
2890! 

2900  !  CALCULATE  AVERAGE  TEMPERATURE  OF  INNER  HEATED  HALL. 

2910! 

2920  Sum-0. 

2930  IF  I  type-0  THEN 
2940  FOR  1-40  TO  69 
2950  Sum-Sum+T(I) 

2960  NEXT  I 
2970  Twi-Sum/30 
2980! 

2990  !  CORRECTION  FOR  BROKEN  T/C  CHANNEL  47  IN  INNER  CURVED  HEATER. 

3000! 

3010  ELSE 

3020  FOR  1-40  TO  69 
3030  IF  1-47  THEN  3050 
3040  Sum-Sum-»T  ( I ) 

3050  NEXT  I 
3060  Tui -Sum/29 
3070  END  IF 

3080  PRINT  USING  'I2X.  "INNER  INSULATION;"" 

3090  PRINT  USING  "14X,""T/C  Number:  71  (Tinsil)  72  (Tinsi2)  73  (Tinsi3)" 

N  ■ 

3100  PRINT  USING  "14X."-Tenp  (C):  .3(3D.DD.8X) " : T ( 70 ) , T( 71 ) ,T ( 7" ) 

3110  Tinsi1-T(70) 


3120 

3130 

3140 

3150 

3160 

3170 

3180 

3190! 

3200 

32101 

3220 

3230 

3240 

3250 

3260 

3270! 

3280 

3290! 

3300 

3310! 

3320 

3330! 

3340 

3350! 

3360 

3370! 

3380 

3390 

3400 

DD," 

3410 

DD," 

3420 

DD,** 

3430 

DD,” 

3440 

DD,” 

3450 

DD,” 

3460! 

3470 

3480! 

3490 

3500 

3510 

3520 

3530 

DD,” 

3540! 

3550 

3560! 

3570 

3580 

DD,” 

3590! 

3600 

3610! 

3620 

3630 


Tinsi2«T(71 1 
Tinsi3-T(72) 

PRINT 

PRINT  USING  ■12X,”0UTLET  TEMPERATURE:*” 

PRINT  USING  ■14X.**T/C  Number:  74  75  76  77*** 

PRINT  USING  ’MX, ••Temp  (C):  •*.4(3D.DD.2X) * ;T( 73) , T(74) , T(75) ,T(7G) 

PRINT  USING  •«,#• 

!  CALCULATE  AVERAGE  OUTLET  TEMPERATURE. 

Sum~0 . 

FOR  1-73  TO  76 
Sum-Sum+T(I) 

NEXT  I 
Tout-Suffl/4 

!  CALCULATE  DIFFERENCE  BETWEEN  CHANNEL  INLET  AND  OUTLET  TEMPERATURES. 
Tdlf f-Tout-Tin 

!  CALCULATE  CHANNEL  BULK  TEMPERATURE. 

Tbllr-(T  in+Tout)/2 
!  PRINT  CALCULATED  TEMPERATURES. 

PRINT  USING  •9X,**THE  FOLLOWING  TEMPERATURES  WERE  CALCULATED: *•’ 

PRINT 

PRINT  USING  *12X, **Average  Outer  Wall  Temperature  (Two)  -  *’, 

(C)**“;Two 

PRINT  USING  *12X,”Average  Inner  Wall  Temperature  (Twi)  •  ”, 

(C)’*“;Tui 

PRINT  USING  ■12X,* "Average  Outlet  Temperature  (Tout)  ■  **, 

(C)*"*;Tout 

PRINT  USING  • 12X, ""Average  Inlet  Temperature  (Tin)  “  *", 

(C)""";Tin 

PRINT  USING  •12X,"*Channel  Inlet  and  Outlet  Temp  Diff  (TdifF)  ■  *", 
(C)""";Tdiff 

PRINT  USING  •12X,"* Average  Bulk  Temperature  (Tblk)  -  ”, 

(C)""*;Tblk 

!  CALCULATE  DELTA  T  (Tdel)  FOR  USE  IN  CALCULATING  H  LATER. 

Tdelo-Two'Tblk 
Tdel i-Tui"Tblk 
IF  Ipos-O  THEN  Tdel-Tdelo 
IF  Ipos-1  THEN  Tdei-Tdeli 

PRINT  USING  •12X,*"Mean  Temperature  Difference  (Tdel)  -  "", 

(C)""";Tdel 

!  CALCULATE  DIFFERENCE  BETWEEN  OUTER  AND  INNER  HALL. 


••,3D. 


•■,3D. 

••.3D. 

”,3D. 


••,3D. 


Tudif f-ABS(Tuo-Twi ) 

PRINT  USING  "12X, ““Outer  and  Inner  Wall  Temp  Difference  (Tudiff)  -  "“,30. 
(C)""";Tudiff 

I  CALCULATE  LOCAL  CHANNEL  TEMPERATURES  AT  INLET.  MIDDLE  AND  OUTLET. 

Tcin“Tin+(  .5'»Tdiff/12) 

Tcmid-T in+( 4.5*Tdif  f /1 2) 


V.  y 


3640  Tcout'T in+ ( 1 1 . 5*Tdi f f / 1 2 ) 

36B0  PRINT  USING  ■12X, •"Average  Local  Channel  Temp  at  Inlet  (Tcin)  -  "■,30. 
DD."‘  (C)""*;Tcin 

3660  PRINT  USING  •12X, ""Average  Local  Channel  Temp  at  Middle  (Tcmid)  -  "",30. 
DD.""  (C)""";Tcmid 

3670  PRINT  USING  •12X,*" Average  Local  Channel  Temp  at  Outlet  (Tcout)  -  "",30. 

DD.""  (C)""";Tcout 

3680! 

3690  I  CALCULATE  AVERAGE  LOCAL  HEATED  WALL  TEMPERATURE  AT  INLET. 

3700! 

3710  Sum-0. 

3720  Iadd-0 

3730  IF  Ipos-1  THEN  Iadd-40 
37401 

3750  !  CORRECTION  FOR  BROKEN  T/C  CHANNELS  5.  AND  47  CURVED  HEATER  ONLY. 

3760! 

3770  IF  Itype-1  THEN 
3780  FOR  I-O  TO  6 

3790  IF  (Ipos-0  AND  1-5)  OR  (Ipos-1  AND  1-6)  THEN  3810 
3800  Sum-Suffl-rTl  I  +  Iadd) 

3810  NEXT  I 
3820  Twin-Sum/6 
3830  ELSE 
3840  FOR  I-O  TO  7 

3850  IF  (Ipos-0  AND  1-7)  OR  (Ipos-1  AND  1-6)  THEN  3870 
3860  Sum-Sum+T( I+Iadd) 

3870  NEXT  I 
3880  Tuin-SuB/7 
3890  END  IF 

3900  PRINT  USING  •12X,""Avg  Local  Heated  Wall  Temp  at  Inlet  (Turn)  •  "".SD. 

DD."'  (C)““*:Tuin 

3910! 

3920  !  CALCULATE  AVERAGE  LOCAL  HEATED  WALL  TEMPERATURE  NEAR  MIDDLE. 

3930! 

3940  Sum-0. 

3950! 

3960  !  CORRECTION  FOR  BROKEN  T/C  CHANNEL  16  CURVED  HEATER  ONLY. 

3970! 

3980  IF  Itype-1  AND  Ipos-0  THEN 
3990  FOR  1-9  TO  15 
-(000  IF  1-10  THEN  4020 
4010  Sum-Sum+T(I) 

4020  NEXT  I 
4030  Twmid-SuiT'/6 
4040  ELSE 
4050  FOR  1-9  TO  16 

4060  IF  (Ipos-0  and  1-10)  OR  (Ipos-1  AND  1-9)  THEN  4080 
4070  Sum-Sum+T( I+Iadd) 

4080  NEXT  I 
4090  Tumid-Suin/7 
4100  END  IF 

4110  PRINT  USING  •12X.""Avg  Local  Heated  Wall  Temp  Near  Middle  (Tumid)  -  **,3D 

DD."*  (C)*"*;Twmid 

4120! 

4130  !  CALCULATE  AVERAGE  LOCAL  HEATED  WALL  TEMPERATURE  AT  OUTLET. 

4140! 

4150  Sum-0. 

4160! 

4170  !  CORRECTION  FOR  BROKEN  T/C  CHANNEL  26  CURVED  HEATER  ONLY. 

4180! 
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4190  IF  Itype-1  AND  Ipos-0  THEN 
4200  FOR  1-23  TO  29 
4210  IF  I-2G  THEN  4230 
4220  Sum-Suin+T(  I) 

4230  NEXT  I 
4240  Twout-Sun/6 
4250  ELSE 

42S0  FOR  1-23  TO  29 
4270  Sum-Sum+Tl I+Iadd) 

4280  NEXT  I 
4290  Tuout-Sum/7 
4300  END  IF 

4310  PRINT  USING  •12X,'’Avg  Local  Heated  Hall  Tei»p  at  Outlet  (Tuout) 

DD."  (C)**“;Tuout 
4320  PRINT 
4330  PRINT 
4340! 

4350  I  ENTER  ORIFICE  PRESSURE  DATA. 

4360! 

4370  BEEP 

4380  INPUT  'ENTER  PATH  (inHg).  DPM(inH20).  PlM(inH20).  REl * .Patm.Dpe.Plm.Rel 
4390  PRINT  USING  •9X,"THE  FOLLOWING  ORIFICE  DATA  HERE  ENTERED;*” 

4400  PRINT  ^  .. 

4410  PRINT  USING  •12X,  "PatmC inHg)  DPM(inH20)  PlM(inH20)  REl'” 

4420  PRINT  USING  '1 4X,2(2D.DD,  10X) ,  (2D.DD,7X)  ,5D.D* ;Patii..Dp«,P1».Re1 
4430  PRINT 
4440! 

4450  !  CONVERT  PRESSURE  READINGS  TO  SI  UNITS  AND  PIN  TO  ABSOLUTE. 

4460! 

4470  Patm-Patra*337G.8 
4480  Pdel-Dpiii*248.84 
4490  Pl-Patra-(Plm*248.84) 

4500  Torf-T(33) 

4510! 

4520  !  PHYSICAL  PROPERTIES  AND  CONSTANTS. 

4530! 

4540  R-28G.987  IFOR  AIR, 

4550  Rho-Pl/(R*(Torf+273.15))  IDENSITY  OF  AIR  AT  ORIF  PLATE. 

4560  Gainma-1.40  !F0R  AIR. 

4570  Cp-1006  !Cp  GOOD  FOR  Tbllc  BETWEEN  12  AND  33  (DEG  C). 

4580  Gc-1.0 

4590  Signia-5.669E-8 

4600  Ecu-. 12  lEMISSIVITY  OF  COPPER. 

4610! 

4620  !  DIMENSIONS  OF  CHANNEL. 

46301 

4640  Dc-. 006350 
4650  Hid-. 2540 
4660  Ri-.297 

4670  Pwet-2*(Dc>Hld) 

4680  Ac-Dc*Wid 
4690  Dhd“4*Ac/Pwet 
4700  Apl-. 0774192 
4710! 

4720  !  PROPERTIES  OF  INSULATION. 

4730! 

4740  Xins-.0127 

4750  Kins-3.8E-2 
4760! 

4770  !  ENTER  ORIFICE  CONFIGURATION. 


-  ”.30, 


.'.V.V, 
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4790  BEEP 

4800  INPUT  "SELECT  DIAMETER  OF  ORIFICE  (0-.5334,  1-1.0755  ( inches) )" ,Sue 
4810  IF  Size-O  THEN  Dorf-.013548 
4820  IF  Size-1  THEN  Dorf-. 027318 
4830! 

4840  !  DIMENSIONS  OF  PIPE  AND  ORIFICE  PLATE. 

4850! 

48G0  A-(PI*Dorf2)/4 
4870  Dpipe-.051772 

4880  Apipe-PI*Dpipe‘2/4 
4890  Beta-Dorf /Dpipe 

4900  PRINT  USING  •12X, ■"Dorf (n)  A(m-2)  Beta"’* 

4910  PRINT  USING  ’12X. (Z.GD.GX) , (Z.3DE,GX) .Z.4D*;Dorf .A.Beta 
4920  PRINT 
4930  PRINT 
4940! 

4950  !  CORRELATION  FOR  EXPANSION  FACTOR  BASED  ON  ID  AND  1 /2D  TAPS. 

49G0! 

4970  Y-1-( .41+.35*Beta*4)*(Pdel/(Ga«i»a*P1 )) 

4980! 

4990  !  CALCULATION  OF  FLOW  COEFFICIENT  (K)  IN  SI  UNITS. 

5000! 

5010  B-. 0002+2. 794E-5/Dpipe+(. 0038+1. 01GE-5/Dpjpe)*(Beta‘2+(1G. 5+1 .9B8504E+2*Dp 

ipe)*Beta* 1G) 

5020! 

5030  I  K1  AND  K2  ARE  DUMMY  VARIABLES  TO  CALCULATE  Ko. 

5040! 

5050  K1-.S014-5.3974064E-3/Dpipe*.25 

5050  K2-(.37G+2.8971138£-2/Dpipe*.25)*<».6129E-7/(Dpipe‘2*Beta'2+G.35E-5*Dpipe) 

+Beta‘4+1 .5*Beta‘ 16) 

5070  KO-K1+K2 

5080  K-Ko+1000*B/Rer .5 

5090! 

5100  !  Mu  CORRELATION  GOOD  FOR  Torf  BETUEEN  17  AND  41  DEGREES  CELSIUS. 

5110! 

5120  Mu-4.G971429E-8*Torf+1 .7194722E-5 

5130! 

5140  !  CALCULATE  Mdot  AND  REpipe  AND  COMPARE  TO  PREDICTED  RE1. 

5150! 

5160  Mdo  t-K+A*  Y* ( 2*Gc*Rho*Pde 1 ) ‘ . 5 
5170  Rep 1 pe- ( Mdo  t*Dp i pe ) / <  Mu* Ap i pe ) 

5180  Diff-(ABS(Re1-Repipe)/Re1 )*100 

5190  Rel-Repipe 

5200  IF  Diff>.001  THEN  5080 

5210  PRINT  USING  "9X.”THE  FOLLOWING  DATA  MERE  CALCULATED:”* 

5220  PRINT 

5230  PRINT  USING  *12X, ’"Orifice  Expansion  Factor  (Y)  •  ”,Z.4D,";Y 

5240  PRINT  USING  •12X, ’"Orifice  Floe  Coefficient  (K)  -  ”  .Z.4D, " JK 

5250  PRINT  USING  ’  12X, ’’Density  Based  on  Torf  (Rho)  -  ”,Z.4D,”  (lcg/i|- 

3)”’;Rho 

5260  PRINT  USING  ’ 12X. ’’ViscosJ ty  Based  on  Torf  fMu)  -  ’’.2Z.3DE.”  (kg/m.s) 

’”;Mu 

5270! 

5280  I  Mu  CORRELATION  GOOD  FOR  Tblk  BETWEEN  17  AND  41  DEGREES  CELSIUS. 

5290! 

5300  Mu-4.6971429E-8*Tbllc  +  1 .7194722E-5 

5310  Red- (Mdo t*Dc ) / (Mu*Ac ) 

5320  Rehd-(Mdot«Dhd)/(Mu*Ac) 

53301 


5340  !  CALCULATION  OF  POWER  INTO  HEATER  PLATE. 

5350 ! 

5360  Iadd-0 

5370  IF  Ipos-1  THEN  Iadd-39 

5380  IF  Ipos-0  THEN  Rpr-2. 01729 

5390  IF  Ipos-1  THEN  Rpr-2. 00839 

5400  Qp- ( E»f ( 38+Iadd ) *E*f ( 39+ ladd ) ) /Rpr 

5410! 

5420  I  CALCULATION  OF  HEAT  INPUT  TO  AIR. 

54301 

5440  Oai r-Mdot*Cp*( Tout-T i n) 

5450! 

5460  !  CALCULATE  LOCAL  HEAT  TRANSFER  COEF.  AND  NUSSELT  NUMBER. 

5470! 

5480  Ka 1 r-7 . 7257 1 43E-5*Tcin+ . 024 1 65836 
5490  Hin-Oair/ ( Apl*( Tuin-Tcin) ) 

5500  Nuin-(Hin*Dhd)/Kair 

551 0  KaiT-7 . 72571 43E-5*Tciild+ .  0241 65836 

5520  Hmid-Oair/C Apl*(Tumid-Tcmid) ) 

5530  Nunid-(Hnxd*Dhd)/Kair 

5540  Kair-7. 72571 43E-5*Tcout+. 024165836 

5550  Hout-Qair/{Apl*(Tuout-Tcout) ) 

5560  Nuout-(Hout*Dhd)/Kair 
5570! 

5580  !  Kaxr  CORRELATION  GOOD  FOR  Tblk  BETWEEN  17  AND  41  DEGREES  CELSIUS. 

5590! 

5600  Kair-7. 7257 1 43E-5*Tbll(  + .  024 1 65836 
5610! 

5620  !  CALCULATE  AVERAGE  HEAT  TRANSFER  COEF.  AND  NUSSELT  NUMBER. 

5630! 

5640  Hawg-0air/( Apl*Tdel ) 

5650  Nuavg-(Havg*Dhd)/Kair 
5660! 

5670  !  CALCULATE  HEAT  LOSSES. 

5680! 

5690  Olol -(Kins*Apl*( T insol -T inso2) )/Xins 

5700  Qlo2-(Kin8*Apl*(Tln8o2-Tin8o3) l/Xins 

5710  01o-(01oI+Qlo2)/2 

5720  Qli1-(Kins*Apl*(Tin8i1-Tinsi2))/Xins 

5730  01i2- (Kins*Apl*(T in8i2-T insi3) )/Xins 

5740  01i-(Qlil+01i2)/2 

5750  Rr-2. 0*( ( 1 -Ecu) /( Apl*Ecu ) )+( 1 .0/Apl ) 

5760  Qr-(Signa*(A3S( (T«o+273.15) •4-(Tui+273. 15) '4) ) )/Rr 

5770  Odel-Qp-Oair-Olol-Qr-Qlil 

5780! 

5790  !  CALCULATE  DEAN  NUMBER. 

5800! 

5810  De-.5*Red*((Dc/2)/Rl)'.5 

5820  PRINT  USING  •12X, ■’Viscosity  Based  on  Tblk  (Mu)  -  ■■.2Z.3DE.*’  (kg/i».s) 
••■;Mu 

5830  PRINT  USING  •12X, ■•Therm  Cond  Based  on  Tblk  (Kair)  -  ’’.Z.SDE,”  (W/m.K)’ 

■“iKair 

5840  PRINT 

5850  PRINT  USING  ■12X."Patm  -  •■.6D.DD,”  (N/m  *2)  ’  ’  ‘  :Pa  tm 

5860  PRINT  USING  ■12X.”Pdel  -  ■■.4D.DD,”  (N/m*2)  ’  *■  ;Pdel 

5870  PRINT  USING  ■12X,”P1  -  ■’.SD.DD,”  (N/m‘2)  ■"’ ;P1 

5880  PRINT  USING  ■12X.”Mdot  -  ■■,Z.4D.”  (kg/s)”“  ;Mdot 

5890  PRINT  USING  ’ 1 2X. ’ ’Repipe  -  ■' .6D.D’ :Repipe 

5900  PRINT  USING  ■12X,’"Red  -  ’".SD.D’lRed 

5910  PRINT  USING  "12X.”Rehd  -  ”  .GD.D"  :Rehd 


5930  !  ONLY  PRINT  DEAN  NUMBER  WHEN  HEATING  CURVED  SECTION. 

5940! 

5950  IF  Itype-1  THEN  PRINT  USING  •12X.’*De  -  ■•,4D.D-:De 

5980  PRINT 

5970  PRINT  USING  •12X,"Qp  -  •■.3D.3D,"‘  (Watts)  *•  .7X,  " ‘Qair  -  ■’,30. 3D.”  (W 

atts)’”  ;Qp,Qair 

5980  PRINT  USING  ■12X,“Qlo1  -  ■•.2D.3D,‘’  (Watts)  ‘  "  .7X.  "  "Ql  i  1  -  •■.2D. 3D,” 

(Watts)’*’,-Qlo1 .0111 

5990  PRINT  USING  •12X,”Qlo2  -  ”,2D.3D.”  (Watts)  ”  .7X.  ”01 12  -  ”.2D.3D.*’ 

(Watts)”‘;01o2.Qll2 

6000  PRINT  USING  ■12X.”Qlo  -  ”.2D.3D.’*  (Watts) ’•,  7X,  ”01  i  -  ”.2D.3D,*‘ 
(Watts)’”;01o.011 

6010  PRINT  USING  ■12X,”Qr  -  •*,20. 3D.**  (Watts) ** .7X. **0del  -  **,2D.3D,** 

(Watts)***;Or,Odel 

6020  PRINT 

6030  PRINT  USING  *12X.**Hin  -  **.3D.DD.**  (W/ii|-2C)  **  .7X,  *  *Nuin  -  **,3D.DD,*;H 

in.Nuin 

6040  PRINT  USING  *12X.**H(ind  -  **,3D.DD.**  (W/ffi*2C) ** . 7X, **Numid  -  **.3D.DD,*;H 
mid.Numid 

6050  PRINT  USING  •12X,**Hout  -  **.3D.DD.**  (W/iii'2C)  *  *  .  7X,  *  "Nuout  -  **,3D.DD.*;H 
out.Nuout 

6060  PRINT  USING  *12X,**Havg  -  **,3D.DD.**  (W/iii-2C)  *  *  ,  7X,  *  *Nuavg  -  •*,3D.DD,*:H 

avg.Nuavg 

6070  PRINT 

6080  PRINT 

6090  BEEP 

6100  INPUT  *WILL  THERE  BE  ANOTHER  RUN?  ( 1 -Y.O-N) * .Go  on 
6110  IF  Go.onOO  THEN 
6120  PRINT  USING  *«,#* 

6130  IF  Im-O  THEN  WAIT  600 
6140  GOTO  1700 
6150  END  IF 

6160  PRINT  USING  *8X.**END  OF  RUN*** 

6170  PRINT  USING  *«,#* 

6180  PRINTER  IS  1 
6190  ASSIGN  «File  TO  * 

6200  END 
6210! 

6220  !  CONVERTS  EMF  TO  DEGREES  CELSIUS. 

6230! 

6240  SUB  Tv8v(V.T) 

6250  COM  /Co/  D(7).Aa(76.2) 

5260  Sum~0 

6270  FOR  I-O  TO  7 

6280  Suiii“Suin+D(  I ) *V' I 

6290  NEXT  I 

6300  T-Suii 

6310  SUBEND 

6320! 

6330  !  CALIBRATES  T/C  READINGS. 

6340! 

6350  DEF  FNTcorrd.I) 

6360  COM  /Co/  D(7),Aa(7G,2) 

6370  Tc-Aa(I.O) 

6380  FOR  J-1  TO  2 
6390  Tc*Tc+Aa ( I , J ) *T ‘ J 
6400  NEXT  J 
6410  RETURN  Tc 
6420  FNEND 


APPENDIX  B:  SAMPLE  CALCULATIONS 


Figure  36  below,  shows  the  major  heat  transfer 
components  for  each  of  the  test  sections.  The  sample 
calculations  that  follow,  demonstrate  the  methods  used  by 
the  computer  to  calculate  these  components  as  well  as  the 
Reynolds  number,  average  heat  transfer  coefficient,  and 
average  Nusselt  number  for  each  set  of  data.  Part  A  of 
Appendix  B  contains  a  sample  printout  of  a  data  run, 
followed  by  a  summation  of  the  data.  The  sample 
calculations  are  for  the  concave  curved  section,  but  those 
for  the  straight  section  are  similar. 


Figure  36.  Energy  Balance  in  Test  Section. 
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SAMPLE  CALCULATION  DATA 

m 

THE  FOLLOWING  CHANNEL  CONFIGURATION  WAS  SELECTED: 

Heating  CURVED 

Test  Section. 

Heating  OUTER  Plate. 

i 

Month , 

Day,  and 

T  ime: 

10:17:13:47: i 1 

THE  FOLLOWING  DATA 

WERE  RECORDED: 

Heater  Voltaoe 

(Vh) 

-  61. 

337  (V) 

Resistor  Voltage  (Vr) 

-  5. 

195  (V) 

Precision  Resistor  (Rpr)  “  2. 

01729  (Ohms) 

ft 

OUTER  PLATE; 

T/C  Number: 

1 

2 

3 

4 

5 

Temp  ( C ) : 

42.75 

41  .57 

44.39 

41  .23 

44.55 

T/C  Number: 

7 

8 

9 

10 

1 1 

Temp  (C) : 

46.54 

42.61 

43.51 

43.47 

44.77 

T/C  Number: 

13 

14 

15 

16 

17 

A 

Temp  (C) : 

43.40 

44.18 

44.10 

41  .97 

38.06 

T/C  Number: 

19 

20 

21 

22 

23 

Temp  (O: 

45.79 

48.45 

49.18 

33.83 

45.32 

T/C  Number: 

25 

26 

27 

28 

29 

Temp  ( C )  : 

45.30 

46.47 

30.31 

45.80 

45.98 

OUTER  INSULATION: 

T/C  Number: 

31  ( T insol ) 

32  (Tinso2)  33 

( T inso3 ) 

Temp  (C) : 

49. 

58 

37.29 

29.36 

ORIFICE  TEMP  (Torf)  - 

24.16  (C) 

INLET  TEMPERATURE: 

\* 

T/C  Number: 

35 

36 

37 

38 

a 

Temp  ( C) : 

20.55 

20.48 

20.47 

20.43 

INNER  PLATE: 

T/C  Number: 

41 

42 

43 

44 

45 

Temp  ( C ) : 

21  .43 

21 .31 

21.16 

21.14 

21 .05 

T/C  Number: 

47 

48 

49 

50 

51 

Temp  ( C ) : 

21.03 

22.39 

21.10 

21.05 

21  .38 

1 

T/C  Number; 

53 

54 

55 

56 

57 

Temp  ( C ) : 

21  .23 

21 .27 

21.17 

21.19 

21.15 

T/C  Number: 

59 

60 

61 

62 

63 

Temp  (C) : 

21  .25 

21 .56 

21.47 

21  .74 

21  .66 

T/C  Number; 

65 

66 

67 

68 

69 

Temp  (C): 

22.05 

21  .91 

21  .85 

21.78 

21.74 

i.'  ^ 

» 

INNER  INSULATION: 

T/C  Number; 

71  (T 

1  ns  1 1) 

72  (Tinsi2)  73 

( Tinsi3) 

V  ' 

Temp  ( C ) : 

21 

.53 

21.49 

21  .32 

'-V- 

OUTLET  TEMPERATURE: 

^  •* 

^V- 

T/C  Number; 

74 

75 

76 

77 

^'C^ 

Temp  (C) : 

24.29 

24.38 

24.38 

24.25 

THE  FOLLOWING  TEMPERftTURES  WERE  CALCULATED: 


Average  Outer  Wall  Temperature  (Two) 

Average  Inner  Wall  Temperature  (Tui) 

Average  Outlet  Temperature  (Tout) 

Average  Inlet  Temperature  (Tin) 

Channel  Inlet  and  Outlet  Temp  Diff  (Tdiff) 
Average  Bulk  Temperature  (Tblk) 

Mean  Temperature  Difference  (Tdel) 

Outer  and  Inner  Wall  Temp  Difference  (Tudiff) 
Average  Local  Channel  Temp  at  Inlet  (Tcin) 
Average  Local  Channel  Temp  at  Middle  (Tcmid) 
Average  Local  Channel  Temp  at  Outlet  (Tcout) 
Avg  Local  Heated  Wall  Temp  at  Inlet  (Twin) 

Avg  Local  Heated  Wall  Temp  Near  Middle  (Tumid) 
Avg  Local  Heated  Wall  Temp  at  Outlet  (Twout) 


44. 7S  (C) 
21.46  (C) 
24.32  (C) 
20.48  (C) 
3.84  (C) 
22.40  (C) 
22.36  (C) 
23.30  (C) 
20.64  (C) 
21.92  (C) 
24.16  (C) 
43. B1  (C) 
43.50  (C) 
46.06  (C) 


THE  FOLLOWING  ORIFICE  DATA  WERE  ENTERED: 


Patm( inHg) 
29.79 


DPM( inH20) 
18.58 


P1M(inH20)  RE1 

12.55  15500.0 


Dorf(m)  A(m*2)  Beta 

0.027318  5.8G1E-04  0.5277 


THE  FOLLOWING  DATA  WERE  CALCULATED: 


Orifice  Expansion  Factor  (Y) 
Orifice  Flou  Coefficient  (K) 
Density  Based  on  Torf  (Rho) 
Viscosity  Based  on  Torf  (Mu) 
Viscosity  Based  on  Tblk  (Mu) 
Therm  Cond  Based  on  Tblk  (Kair) 


0.9852 

0.6323 

1.1424  (kg/m*3) 

18.330E-06  (kg/m.s) 
18.247E-06  (kg/m.s) 
2.590E-02  (W/m.K) 


Patm 

100594.87 

(N/m‘2) 

Pdel 

4623.45 

(N/m-2) 

PI 

97 '*71 .93 

(N/m*2) 

Mdot 

U.0375 

(kg/s) 

Repipe 

50352.2 

Red 

8097.1 

Rehd 

15799.1 

De 

418.6 

Op 

m 

157.971 

(Watts) 

Qair  • 

144.950 

(Watts) 

Qlol 

M 

2.847 

(Watts) 

Qlil  - 

.009 

(Watts) 

Qlo2 

m 

1 .838 

(Watts) 

Qli2  - 

.038 

(Watts) 

Qlo 

m 

2.342 

(Watts) 

Qli  - 

.023 

(Watts) 

Qr 

» 

.751 

(Watts) 

Qdel  “ 

9.415 

(Watts) 

Hin 

« 

81  .90 

(W/m‘2C) 

Nuin 

-  39.39 

Hmid 

m 

86.78 

(W/m'2C) 

Numid 

-  41.58 

Hout 

m 

85.51 

(W/m'2C) 

Nuout 

-  40.70 

Havg 

m 

83.75 

(W/m'2C) 

Nuavg 

-  40.07 

H 

s 

61.337  V 

PR 

= 

5.195  V 

PR 

= 

2.0173  a 

atm 

= 

29.79  in.  Hg 

P 

18.58  in.  H2O 

1 

= 

12.55  in.  H2O 

in 

= 

20.48  C 

out 

= 

24.32  C 

wo 

= 

44.76  C 

wi 

= 

21.49  C 

ins , i  1 

= 

21.53  C 

ins ,  i  2 

= 

21.49  C 

ins ,  01 

= 

49.58  C 

ins , 02 

= 

37.29  C 

orf 

=: 

24.16 

air 

= 

25.90  X  10“^  W/mC 

pair 

= 

1006  J/Kg  K 

air 

= 

4.6971  X  10“®  X  T^.  +  1.71947  x  10 

ins 

3. 

80  X  10"^  W/mC 

^ins 

= 

0.0127  m 

cu 

= 

0.12 

5.669  X  10~®  W/ra^K^ 

= 

.5277 

= 

1.  40 

= 

1  Kg  m/Ns^ 

wo-wi 


^pipe 


D  . 
pipe 


=  286.98  Nm/Kg  K 


=  1.0 


=  .07742 


=  .0016 


=  .00211  m' 


=  .00635  m 
=  .0518  m 


=  .027318  m 


=  .5207  m 
=  .297  m 


2.  TEMPERATURE  CALCULATIONS 


a.  Bulk  Temperature 


T  '"out  ,  20^48  .  24,32  .  22.40  C 

Dl  K  2  2 


b.  Mean  Temperature  Difference  (AT) 


AT  =  T_  -  T.  ..  =  44.76  -  22.40  =  22.  36 

wo  bl  k 


c.  Temperature  Difference  Between  Outer  and  Inner  Wall 
^"^wdiff  ^ 


=  T.  -  T  .  =  44.76  -  21.49  =  23.  27  C 
wdiff  wo  wi 


REYNOLDS  NUMBER  CALCULATIONS 


a.  Re„,„  -  - (.03751  (.0516)  . 

''plpe'alr  (-00211)  (18. 33x10-^ 


m  D 


b.  Re^  = 


A  u  . 
c  ^ai r 


(.0375)  (.00635)  ^  8097.1 

(.0016)  (18. 25x10“®) 


m  D 


c .  Re 


hd 


hd  Ay. 

c  air 


4  X  A 


(1)  D 


hd  P, 


wet 


c  _  (4)  (.0016)  _ 

- 752^7 - .01229  m 


m  D 


Re 


hd 


hd 


A  |i  . 
c^ai  r 


^  (.0375)  (.01229)  ^  15799.1 

(.0016) (18. 25x10“®) 


HEAT  CONVECTED  TO  AIR  CALCULATION 

a.  Q,{^=  mC„(T„.  .-T.  )  =  (. 0375)  (1006)  ( 24. 32-20. 48) =144 . 95  W 
oiiT  p  olic  in 


AVERAGE  HEAT  TRANSFER  COEFFICIENT  CALCULATION 


h  = 


^ai  r 


144.95 


Apj^  AT  (.07742)  (22.  36) 


=  83.72  W/m  C 


AVERAGE  NUSSELT  NUMBER  CALCULATION 

—  _  ^  °hd  _  (83.73) (.01229)  _ 

ha  -  KITT  -  25.90  X  10-3 


39.73 


/  f 


APPENDIX  C:  EXPERIMENTAL  UNCERTAINTY 


The  uncertainties  for  the  major  variable  in  the 
experiments,  were  calculated  in  accordance  with  the  method 
described  by  S.  Kline  and  F.  McClintoch  [Ref.  33].  The 
estimates  of  the  uncertainty  in  the  measured  quantities  were 
made  conservatively.  As  a  result,  there  is  considerable 
confidence  in  the  uncertainties  as  calculated. 

The  following  equations  were  used  to  calculate  the 
uncertainties: 


A.  Aa  2  dP,  2  dT^^p  2 

(1)  ((f^)  +  (p-^)  +  (t~^)  ) 

1  orf 


dm  _  r,,dY,^  .  .dK.^  ^  /dA.^. 

(2)  -  =  [  (  (^)  +  (j^)  +  (j^)  ) 


dp  •  2  Had  2 

+  (1/4  +  (fl^)  )] 

^ai  r 


dQ^.  dm  2  d  2  d(T^,.-T.  )  2 

(3)  =  ((— )  +  (-r  ^  )  ) 


Qair 


C 

pair 


T  -T . 

out  in 


aT  dQ  .  2  dAp.  2  2 

(4)  +  (r-^)  +  (^)  ) 

h  ^air  ''PL 


dNu.  ,  -  _  dD.  ,  2  dK^.  2 

(5)  =  ((^)^  +  (g-^)  +  (jr^)  ) 

Nu^^  h  hd  air 


<  V-  •  • ' 

v'  'X*- 


The  following  quantities  had  their  uncertainty 


calculated  by  dividing  the  estimated  error  in  the  quantity 
by  the  value  of  the  quantity,  and  these  uncertainties  are 
considered  constant  for  the  range  of  values  in  this 
experiment ; 


Quanti  t 


Uncertaint 


Quantity 

Error 

"^blk 

+  .  1°C 

'^in 

+  .  1°C 

^orf 

+  .  1°C 

"^out 

+  .  1°C 

T 

wo 

+  .1°C 

T  -T . 

out  in 

+  .  2°C 

AT 

+  .2°C 

For  hydraulic  Reynolds  numbers  greater  than  17000  the 
error  in  the  outlet  temperature  was  much  higher,  as 
mentioned  previously,  and  the  following  errors  apply: 


Quantity 

Error 

T 

+  .  2°C 

blk 

Tj 

+  .  1°C 

in 

T  c 

+  .  1°C 

orf 

T  ^ 

+ 

• 

o 

o 

out 

— 

T 

+  .  1°C 

wo 

T  -T . 

out  in 

+  .4°C 

AT  +.3°C 


The  uncertainties  for  the  following  quantities  were 
calculated  using  the  sample  data  for  the  outer  curved  test 
section  at  a  hydraulic  Reynolds  number  of  15800: 


Quantit 


Uncerta i nt 


The  uncertainty  for  the  following  quantities  were 
calculated  for  the  outer  curved  test  section  at  a  hydraulic 
Reynolds  number  of  21200: 


Quantity 


Uncertainty 


By  multiplying  the  quantity  by  its  uncertainty  the 

possible  error  can  be  calculated.  The  previous  uncertainty 

calculations  for  a  hydraulic  Reynolds  number  of  15800 

yielded  Nu^^  =  40.07  +2.73  and  a  Rsj^^  =  15800  +  928.  For 

the  run  at  a  hydraulic  Reynolds  number  of  21200  the  Nu^^^  = 

44.71  +  5.03  and  a  Re.  ,  =  21200  +  1232. 

hd  — 

The  major  source  of  uncertainty  in  the  average  Nusselt 
number  is  the  uncertainty  in  the  difference  between  the 
outlet  and  inlet  air  temperature.  If  this  uncertainty  can 
be  brought  below  four  percent,  then  the  uncertainty  of  the 
average  Nusselt  number  will  be  dictated  by  the  uncertainty 
of  the  measured  dimensions  of  the  channel.  In  particular, 
the  channel  height  (D^)  which  has  an  uncertainty  = 

.0400.  The  uncertainty  of  the  channel  height  can  only  be 
reduced  by  increasing  the  channel  height,  and/or  using 
materials  with  better  tolerances. 
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